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1 Introduction 

In the recent series of papers [5, 6, 7, 9, 12, 13, 14, 15, 16] a special class of Inonu-Wigner 
contractions was introduced, namely "analytic contractions". For an analytic contraction the 
contraction parameter R, which is the radius of an re-dimensional sphere or pseudosphere Uq + 
eu 2 = R 2 , (e = ±1) figures in the separated coordinate systems, and in the eigenfunctions and 
the eigenvalues for the Laplace-Beltrami operator (or Helmholtz equation). With the help of 
analytic contractions we have established the connection between the procedure of separation 
of variables for homogeneous spaces with constant (positive or negative) curvature and flat 
spaces. For instance it has been indicated how the systems of coordinates on the sphere 52 and 
hyperboloid H2 transform to four systems of coordinates on Euclidean space E 2 . 

The goal of this note is to establish the contraction limit R — > 00 for eigenfunctions (or 
basis functions) of the two-dimensional Helmholtz equation on the two-sheeted hyperboloid 
H 2 : Uq — u\ — u 2 = R 2 , 

Alb* = -^±^*, a = -1/2 + ip, (1) 

where the Laplace-Beltrami operator in the curvilinear coordinates (£ , £ 2 ) has the form 
Id d 

A LB = -~^^7iV99 tk ^k, 9 = I det{g ik )\, g ik g kfl = 6?, 
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with the following relation between local metric tensor gik(£,), {i,k = 1, 2) and the ambient space 
metric G^ v = diag(l, —1, —1), 

9ik(£) = Gfiu d p g t k , i,k = l,2, n,v = 0,1,2. 

It is well known that the Helmholtz equation (1) on the two-sheeted hyperboloid admits sep- 
aration of variables in nine orthogonal systems of coordinates [10, 18]. These nine systems of 
coordinates can be separated into three classes. The first class includes three systems of coordi- 
nates which are of subgroup type: viz. pseudo-spherical, equidistant and horicyclic coordinates, 
the second class includes three non subgroup type of coordinates: semi circular parabolic, elliptic 
parabolic and hyperbolic parabolic. The last two coordinate systems in the general case contain 
the dimensionless parameter 7 and in the limiting cases when 7 — > or 7 — > 00 transform to 
horicyclic or equidistant systems of coordinates. These systems of coordinates have an important 
property; the Helmholtz equation (1) admits the exact solution (which we call exactly solvable) 
in terms of "classical" special functions, namely Bessel and Legendre functions or hypergeo- 
metric one. The third class of coordinates consists of elliptic, hyperbolic and semi hyperbolic 
coordinates each of which belong to the class of non subgroup type. Two of them, elliptic and 
hyperbolic coordinates, also contain a dimensionless parameter which is included in the metric 
tensor and Laplace-Beltrami operator. The Helmholtz equation (1) after the separation of vari- 
ables in these coordinates leads to the two Heun type differential equations with a four singular 
points (which we call non-exactly solvable equations) and whose solutions appear as Lame or 
Lame-Wangerin functions [3, 4, 8]. In contrast to the classical special functions which can be 
appear in terms of hypergeometric functions, the Lame or Lame-Wangerin functions are defined 
in terms of infinite series where the expansion coefficients can not be written in explicit form 
because they are the subject the third or higher order recurrence relations. 

The contraction limit of a basis function is not a trivial task in the case of exactly or non- 
exactly solvable equations. Some calculations come from the papers [5, 6, 7], but many of them 
are not known till now. 

In this paper we restrict ourselves to the contraction limit R — > 00 for four kinds of orthogonal 
basis functions: horocyclic, semi circular parabolic, elliptic parabolic and hyperbolic parabolic, 
which is new. We shall present normalizable eigenfunctions but do not give their normalization 
constants explicitly, except in the case of horocyclic coordinates which is particularly simple. 
Here we have also included for completeness the contraction for pseudo spherical and equidistant 
wave functions (see [7, 12]). 

We hope that our results beside possible application to the theory of special functions will be 
useful also in the investigation of super integrable systems which admit separation of variables 
on the two-dimensional hyperboloid. We think that they can be generalized for the three and 
higher dimensional hyperbolic space when these six systems of coordinates are the sub systems 
of more complicated systems of coordinates and where the procedure of variable separation leads 
to similar differential equations. 

2 Contractions 

2.1 Pseudo-spherical basis to polar 

The first coordinate system is the pseudo-spherical system (r > 0, 92 6 [0, 27r)) 
uo = -Rcoshr, u\ = -Rsinhr cos 99, 112 = i?sinhrsin<^, 
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and the orthogonal basis functions of equation (1) are [7, 8, 12, 13] 

^ w (r^) = ^H 1/2 (coshr)e^ ; 

where Py{z) is the Legendre function. 
In the contraction limit R —> oo we have 

r 

sinhT ~t~— , if — > ip, p~ kR, 
R 

where (r, ip) are the polar coordinates in Euclidean plane. There is now the matter of how this 
contraction affects the basic eigenfunctions that can be computed on the hyperboloid. Using the 
well-known representations of the Legendre function in terms of the hypergeometric function [1], 

M ^ ^ _ r(l/2 + i, + M) / , n N , n M 1 



P |m| fcosh r) - L A±[^P+M} ( sinh n m ( cosh 1 
^P-i/^ coshT >- Y{l/2 + ip-\m\) V mh 2) l C ° Sh 2/ \my. 

x 2-P1 ^1/2 + ip + |m|, 1/2 — ip + \m\; 1 + |m|; — sinh 2 y . 
asymptotic formulae for gamma- functions at large z [1], 



T(z + a) 



r(z + /3) 

and taking into account that at R — >• oo 



11 r 

( — + |m| + ip, — + |m| — ip;l + \m\; — sinh 2 — 



^,2 r 2\ / o \ H 

F 1 I 1 + H; - 



/ 2 \ M 



where J u (z) is the Bessel function [2], we get 

/ 2 \ 



<J 1/2 (c shr) « PS_ 1/2 (^1 + — 2 j « (-fcfl) |m| J ]m \(kr). 

Finally, the pseudo-spherical functions in the contraction limit R — > oo take the form 

*% n (T t <p)ni-kR)Wj H (kr)j m <<>, 

i.e., the pseudo-spherical basis up to the constant factor contracts into polar one. The correct 
correspondence to give limiting orthogonality relations in the polar coordinates r and (p can be 
obtained using the contraction limit of the normalization constant (see [7] ) and the results 

sinh rdrdip — > -—^rdrdip, 5(p — p) — > ~^S(k — k'), 
R R 



and 



f°° i 1 

J J\ m \{kr)J\ m \{k'r)rdr = -6(k - k). 
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(2) 



2.2 Equidistant basis to Cartesian 

The coordinate system is the following one {t\ , T2 € M) 

no = -Rcoshri coshr2, u\ = i?coshri sinhr2, 112 = -Rsinhri, 

From the above definition we have that 

sinhri = — , tanhr2 = — 
it uq 

and in the limit of R — > 00 we get 

sinh n ~ n ~ — , sinh T2 ~ T2 ~ — , 
R R 

where x and y are the Cartesian coordinates in the Euclidean plane arising from contraction 

R — > 00. The orthogonal wave function takes the form [7, 8, 12, 13] 

^ Q (ri,r 2 ) = (coshrO-^Pi^^^-etanhrOexp^rs), 

where e = ±1. 

To perform the contraction we use the formula [1] 

~ 2 F l + v), \{l - u + u); |; z 2 ) 

_r(i(i-(,i + ,))r(i + i(,-^) 
2 FiQ(i-At-i/),i + 5( l/ -M);i;« 2 



v 2\~M/2 



2z- 



r(i(i + i/- M ))r(-i(i/ + ^)) 

and rewrite the Legendre function in terms of the hypergeometric function 



R 1 /n(— etannri) = — j» r 7^ r- 

v r(f-a)r(f-&) 



2^1 



1 1 , 1 , 2 

- — a, — — 0; -; tanh n 



+ 2e tanh n — rf 



r(|-a)r(|-6) p /3 



r(i-a)r(i^y^U" a 'i"' ; 2 ;tanh Tl 

where a = z(p + ^)/2, 6 = i(p — v)/2. In the contraction limit i? — > 00 we put 
Then using also (2) we have the asymptotic formulae 



lim 2F1 

R— >oo 



lim 2F1 

R— >oo 



111 

- - a, - - b; -; tanh 2 n ) = o^i 



- - a, - - b; -; tanh ri ) = o^i 



1 y 2 fc 
2' 4 



2J.2 
2 



COS /C22/, 



3. y 2 kj\ _ 1 



.2' 4 y fc 2 y 
where /c 2 + k\ = k 2 give us the contraction limit for P\^_-^ / 2 (~ etanhri), 



sin /c 2 y, 



jikR 



V 



^^(-etanhn) » ^_i/2 
and finally for vpp;/ 



7T 2 tkR exp(ie/c2y) 



3 _ iR(k+k 1 )\ F / 3 i-R(fc-fci) 



7T2 



3 _ igCfc+fci) \ r / 3 _ ig(fc-fci) 



exp(i/cia; + iek 2 y). 



(3) 



(4) 
(5) 



The last result up to the constant factor coincides with the contraction limit of equidistant wave 
function to Cartesian one in Euclidean space. 
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2.3 Horocyclic basis to Cartesian in E 2 



In the notation of the article mentioned horocyclic coordinates on the hyperboloid can be written 

(x G K, y > 0) 



uq = R 



x 2 + f + l 

2y 



u\ = R 



x 2 + y 2 — 1 
2£ : 



X 

u 2 = R-. 

y 



From these relations we see that 
U 2 R 



X 



y 



Uq — U\ Uq — U\ 

In the limit as R — > oo we obtain 



x 



y_ 
IV 



v -> i + 



x 

R' 



where x and y are the Cartesian coordinates in the Euclidean plane. The horicyclic basis 
functions satisfying the orthonormality condition 



R' 



/OO PC 
dx 
-oo JO 



have the form 



^?{x,y) = f-^^V-yK l Myy s ^ 

where K v {x) is the Macdonald function [2]. 

To effect the correct contraction we further require that p —> kR, and s 
k = \Jk\ + k\. Consequently we need the asymptotic formula for 

K ikR (k 2 (R + x)) 

as R — > oo. For this we use the asymptotic formula [2] 



k<iR where 



K iv {x) 



(y2 _ 3.2)1/4 



exp 



TTV\ 



sin 



( — — \/ v 2 — x 2 + z/cosh 1 — ) , 
\4 xJ 



(6) 



valid if v > x 3> 1 and both v and x are large and positive. If we use this formula for the 
Macdonald functions and take into account that 



kR kirR 

2 6 2 



\J p sinh irp ~ 
then we finally get 




1 k 

■®f°(y,x) ~ —^j — sm(k l x-M)exp(ik 2 y), 



M=^- + R 
4 



cosh 1 /ci 

k 2 



The correct correspondence to give the limiting result in the Cartesian coordinates x and y can 
be obtained using the delta-function contractions 



5(p-p>)5(s-s') 



k 

k x R 2 



<S(fci - k[)5(k 2 - k' 2 ). 
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2.4 Semi circular parabolic to Cartesian coordinates 

These coordinates are given by the formulae (rj, £ > 0) 

uo = R — ^ — ' U1 = R — ^ — ' U2 = R ^r- 

In terms of these coordinates we see that 

2 JR 2 + u\ + u 2 a2 JR 2 + ul - u 2 

V = , £ = • 

Uq — U\ Uq — U\ 

In the limit as R — > oo we have 

9 x + y n x — v 

11 1 + _ fT' ^ 1+ r ' 

We see that in the contraction limit the Cartesian variables x and y are mixed up. To have 
a correct limit we can introduce the equivalent semi-circular parabolic system of coordinate 
connected with previous one by the rotation about axis uq through the angle 7r/4, 

Uq = R — , U\ = it—— — — h 1 



U2 = R 



8£ry 2 V 2£r) 8£r/ 



2 V 2£?7 8£r? J 

In terms of a new coordinates we have that 

2 = y / 2^ 2 + (n 2 -n 1 ) 2 -K-m) ^ 2 = ^2R 2 + {u 2 -u l ) 2 + {u 2 -u l ) 
V2u - (u 2 + ui) V2u - (u 2 + m) 

and in the contraction limit R->oowe obtain 

The suitable set of basis functions are [4, 8] 

for s > and 

*pf P (e,^) = v^7^p (V=if) Jip {V^~s V ) , 

for s < 0. The correct limit is then obtained by choosing s = R 2 (k 2 — k 2 ) and p = kR, where 
k\ + k\ = k 2 . To find the contraction limit of the basis function let us use the asymptotic 
relation for the Bessel function of pure imaginary index [17] 

2TrJ ip (z) ~ — ^ )1/4 exp (i^p 2 + z 2 - ipsinh- 1 V - - i^j exp (^) , 



and for the MacDonald function the formula (6). Then 
Ji P (v^e) ~ JikR (R\/k 2 -k 2 ^J l + V2^j 
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C2 



kR 



2*y/2irk 2 R 



exp ik 2 y + iR 



y2k2 — k sinh 



-i 



7T 



and 



\/2A; cosh 



A; 



Using the last asymptotic formulae it is easily to get for a large R 

^p? P (£, ??) ~ Ry ^^ 2 ex P ( ifc 22/ + **i ~ *f ) sin ( fc i x ~ \ + ' 



where 



<Ji + 5 2 = \[2R(k\ + fca) - flfcsinlr 1 M , 

5i - S 2 = V2R(k 2 -h)- Rk sinh 



-i / fc 2 - fci 

V^2 + fcl. 

These expressions are arrived at under the assumption that k\> k\. In case of k\ < k\ we can 
make the interchanges x with y, and k\ with k 2 . 

2.5 Elliptic parabolic basis to parabolic 

Elliptic parabolic basis contracts into a parabolic one on E 2 . In these coordinates the points on 
the hyperbola are given by [9 S (— 7r/2, 7t/2), a > 0], 



u\ = R 



sinh 2 a — sin 2 
2 cosh a cos 



U2 = -Rtanha tan( 



cosh 2 a 



u Q + ^Jul~ R 2 



_ ^cosh 2 a + cos 2 9 
2 cosh a cos # ' 

From these relations we see that 

co^9 = U °-^- ^, 
In the limit as R — > oo we obtain 

cos 2 9 1 — — , cosh 2 a — > 1 + — , 

it it 
where the parabolic coordinates (£, 7?) are given by 

y = €v, £,v> o. 

The elliptic parabolic wave functions on the hyperboloid have the form [4] 

^!(a,9) = v^i£. 1/2 (suie)J^L 1/2 (tanlia). 
To effect the contraction we take s — > kR and p — > kR and 



tanha 



R 



smtf ~ 



R 



8 



E. Kalnins, G.S. Pogosyan and A. Yakhno 



The separation equation 



da 2 



P 2 + l/4 \ F{a) 
cosh 2 a J 



0, 



becomes 



+ A + k 2 e F(0 = 



where A is now the parabolic separation constant and we impose the condition 



A = R(k 2 — k 2 ) 



or 



k = k + 



2kR 



+ 0(R 



-2\ 



This requires taking the limit of a and the dependent part of the eigenfunctions P?^_ 1 y 2 (tanh a) 

and Pf—\ /2( sm ^)- ^he nm it can be established from the known representation of the Legendre 
function (3). This results in the asymptotic formula 



ikR-h 



-)ikR+ 

t-i / 3 iA 

1 U 4k 



ikR) 



Dir.! 



where D u (z) is a parabolic cylinder function [1] 



D v {z) = 2^ 2 



7re * 



1 



r(i=H) 



i^i 



2' 2' 2 



,(V=2ifc£), 



:V2 



l-v_3_z 2 
2 ' 2' T 



If we look at the theta-dependent part of the eigenfunctions and the corresponding limit taking 
sin = we then obtain the limit 



z>ikR 
ikR+ 



and finally 



r( 



3 _ iA 

4 4fc 



D i 



->2ikR+i 



From these expressions we see that we do indeed obtain the correct asymptotic limit. 



2.6 Hyperbolic parabolic to a Cartesian basis 

Hyperbolic parabolic basis contracts into a Cartesian one on E2. In these coordinates the points 
on the hyperbola are given by [9 G (0, 7r), 6 > 0], 



cosh 2 6 + cos 2 6 

u ° = R n ■ v. 

zsmna smO 

From these relations we see that 



sum b — sin _ . , , 

u\ = R — , U3 = R cot coth 0. 

zsmna smtf 



cos 2 



u - v 7- "? + R2 

Uq - Ui 



cosh 2 b 



up + \A 2 + R 2 

UQ-U\ 
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In the limit as R — > oo we can choose 

,2 



cos 2 6 



y 



cosh 2 6^2(1 + 



2R2 , . ^ , R 

The hyperbolic parabolic basis function on hyperbolid can be chosen in the form [4] 



HP, 



(sinh b sin ^) 1 / 2 P^_ 1/2 (cosh &)P*f_ 1/2 (cos 



1/2 



(7) 



To proceed further with this limit we take p 2 ~ k 2 R 2 and s 2 ~ (k\ — k^)R 2 (the case of s 2 < 0, 
or k\ < corresponds to the discrete spectrum of constant s, and we do not consider this case 
here) where k 2 + k% = k 2 , then using the relation between Legendre function and hypergeometric 
functions (3) and formulae (4), and (5), we obtain for 9 depending part of basis function 



cos 6) ~ P 



ikR 



V 



y/j^R-1/2 \y/2R, 

2 ikR ^exp(ik 2 y) 



iR_ 

2 



k ~\~ *\J k^ k^ 



\/k 2 — k 2 



For the limit of the b dependent part of the eigenfunctions we must proceed differently. In fact 
we need to calculate the limit of 



pikR 
i^/k\-k\R-\l2 



2 1 + 



R 



as R — > oo. We know that the leading terms of this expansion have the form 

Aexp(ikix) + B exp(— ik±x), 
and we now make use of this fact. By this we mean that 



lim P ik %^ 



2|t + 5 



AeyLp(ikix) + B exp(— ik\x) 



where the constants A and B depend on R. It remains to determine A and B. To do this let us 
consider x = 0. We then need to determine the following limit 



lim p ik R ., 



:(V2) =A + B. 

From the integral representation formula [1] 
Y{-u - n)Y{\ + v - ii) 



(8) 



r(i/2-/i) 



i 



-M/2 



(z + coshi 



cosh + 1/2] t) dt, 



(9) 



the above limit requires us to calculate as R — > oo 

/•oo 

/ (\/2 + coshi) 
J o 



AkR-l/2 



RJkf 



2 - k 2 



t dt. 



It can be done using the method of stationary phase [11]. We obtain 



pikR 
-\+iR^Jk\-k\ 



(V2) 



hf(V^^) r{ i_ ikR) 



RUkl-kl + k) T \ + iRUk 2 -k 2 -k 
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(10) 



"1 



By considering the expression for the derivatives of the Legendre function (9) at x = 0, we derive 
the expression 



I™ 



x=0 



then 

Comparing the above relation with (8), we obtain that A = and B is equal to (10), that is 
(cosh 6) ^Be~ ixkl . 

is 2 

Finally, solution (7) contracts as follows 

2 ikR y/TrB 

^ s(6 ' 0) 7 Zl 7 1 f^2\ exp(i/c 2 y - ik lX ) 



3 Conclusion 

In this note we have constructed the contraction limit R — )■ oo for the unnormalized wave func- 
tions which are the solution of the Helmholtz equation on the two dimensional two sheeted 
hyperboloid in four coordinates systems, namely, horocyclic, semi circular parabolic, elliptic 
parabolic and hyperbolic parabolic. Of course the complete analysis of the contraction prob- 
lem must include also the solutions of Helmholtz equation in the three additional systems of 
coordinates as elliptic, hyperbolic and semi hyperbolic. We will study this in the near future. 

We have not presented limits associated with nonsubgroup coordinates. The extension of 
these ideas to problems in higher dimensions is natural and will be presented in subsequent 
work. 
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